Orientation Dependence of Step Stiffness: Failure of SOS and Ising Models to 

Describe Experimental Data 
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We have investigated the step stiffness on Cu(00f ) surfaces as a function of step orientation by 
two independent methods at several temperatures near 300 K. Both sets of data agree well and show 
a substantial dependence of the stiffness on the angle of orientation. With the exception of steps 
oriented along (110), the experimental stiffness is significantly larger than the stiffness calculated 
within the solid-on-solid (SOS) model and the Ising-model, even if next nearest-neighbor interactions 
are taken into account. Our results have considerable consequences for the understanding and for the 
theoretical modeling of equilibrium and growth phenomena, such as step meandering instabilities. 
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Self-assembly of nanostructures is a key route on the 
roadmap of device technology at ever-shrinking lengths. 
Several methods for self-assembly of nanostructures via 
kinetics, kinetic instabilities, or as equilibrium phenom- 
ena are associated with the properties of steps on solid 
surfaces [1]. In the context of nanostructures the step 
stiffness (3 — which describes the resistance of a step to 
meandering — is one of the key parameters in the widely 
applicable step continuum model [2]. As the free energy 
cost per length (along the macroscopic edge direction) 
for elongating the microscopic length of a step due to 
fluctuations, (3 plays a central role in virtually all prop- 
erties involving step excitations. For steps in equilibrium, 
the experimental and theoretical work focussed on spatial 
and temporal fluctuations and the repulsive or attractive 
interactions between fluctuating steps [2, 3]. When steps 
represent a growth front, stable step- flow growth occurs 
due to the presence of an Ehrlich-Schwoebel (ES) barrier 
[4, 5] which hinders interlayer diffusion. This barrier, 
however, can also give rise to a step meandering insta- 
bility: Bales and Zangwill [6] showed that this form of 
instability selects a fastest-growing mode of meandering, 
which determines the morphology of the surface in the 
long-time limit. More recently it was shown that mean- 
dering can also be caused by a kink Ehrlich-Schwoebel 
barrier (KESE) [7, 8]. Experimentally, a meandering in- 
stability was observed on Cu(001) vicinals by Maroutian 
et al. [9]. The existing theories, while accounting qual- 
itatively for the effect, fail to give such quantitative as- 
pects as the temperature-dependence of the dominant 
wave length [9] . 

Regardless of the precise cause of the instability, the 
forces driving the step into meandering compete against 
diffusion currents tending to smooth out gradients in the 
chemical potential \x along steps. The chemical potential 
is the product of (3 and the curvature k. Thus, it is the 
stiffness which keeps deviations from the mean direction 



low. For steps in equilibrium, such as on vicinal surfaces 
in the absence of growth, the stiffness determines the 
magnitude of the equilibrium fluctuations of steps and 
thereby the fluctuations in the terrace- width distribution. 

Unlike the surface stiffness, the step stiffness has a de- 
fined value for any orientation since a step is thermo- 
dynamically rough above temperature T = 0. Present 
analytical theories of the meandering instabilities merely 
consider an angle independent ("isotropic") stiffness [8- 
13] and estimate (3(9) as its value for a step orientation 
along the direction 9 = oi close packing (cf. Eq. (3) be- 
low). In Monte-Carlo simulations (e.g. Ref. [8]), where 
the modeling of the step meandering is based on the solid- 
on-solid (SOS) model [14], a particular angular depen- 
dence of the stiffness is automatically built in by consid- 
ering an angle dependent line tension (step free energy 
per length) (3(9) (since (3(9) = (3(9) + (3" (9)). 

In this Letter we show that not only the standard 
assumption of isotropic stiffness but even the angular- 
dependent step stiffness of the SOS or Ising model [15] is 
far from reality for the much-studied [2, 3] Cu(001) sur- 
face, typical of late-transition/noble metal surfaces and 
on which the meandering instability was observed [9]: 
Our experimentally determined stiffness varies by sev- 
eral orders of magnitude for orientations close to (110), in 
agreement with the SOS and Ising models. For larger an- 
gles, however, the experimental stiffness is about quadru- 
ple the value calculated in the two models [16, 17]. As 
we will argue, this deviation cannot be accounted for by 
including next-nearest-neighbor interactions. Rather it 
must be attributed to many-body interactions such as 
kink-kink interactions and/or effective corner energies. 

Our step-stiffness data come from two independent 
experiments. One set of data is obtained by measur- 
ing the step-step distance correlation function for steps 
which are misoriented by an azimuthal angle 9 from the 
high-symmetry, close-packed (110) direction and there- 



fore contain forced kinks. The second set is obtained 
by analyzing equilibrium shapes of two-dimensional is- 
lands [18, 19]. Here 6 denotes the angle between the 
local island-edge orientation and the (110) direction (so 
that tan 6 corresponds to the local slope) rather than (as 
in Ref. [18], e.g.) the angle denoting the position on the 
perimeter [20]. 

The experiments were performed in a standard ultra- 
high vacuum (UHV) chamber equipped with a variable- 
temperature microscope [21, 22] with high thermal drift 
stability. For the experiments on the step-step correla- 
tion function we used copper surfaces with Miller indices 
(5,8,90). Hence, the surfaces arc tilted relative to the 
(001) surface by the polar angle <f>= 5.98°. The result- 
ing monatomic steps on the surface are rotated with re- 
spect to the close-packed (110) direction by a nominal 
angle 6 = 12.99°, which corresponds to a density 0.23 
of geometric kinks. The accuracy of orientation is natu- 
rally limited by the mosaic structure of the single crystal. 
Crystal cleaning procedures followed established recipes 
[19] after which the sample surfaces revealed clean paral- 
lel steps, which, however, merely on the average display 
the nominal rotation angle of 12.99° of the (5,8,90) sur- 
face. The misalignment 9 of the steps with respect to 
the (llO)-direction varied locally between 0° and 45°. 
The size of the local area of constant angle was of order 
10 4 ran 2 . Fig.l shows an STM image of the Cu(5,8,90) 
surface at 301 K with = 7.13°, i.e. with an average con- 
centration 0.13 of forced kinks. The mean value of 9 of 
the steps was determined by averaging over all steps in 
the image. As is obvious from Fig. 1, the (110) direction 
(black solid arrow in Fig. 1) is easily defined in an STM 
image by the straight step segments between individual, 
monatomic kinks. 

Spatial fluctuations of steps can be analyzed via a step 
correlation function G(y)[25] 

G(y) = ([x(y + y')-x(y')f) = 1 ^\y\a h (1) 

pa\\ 

where y and y' are coordinates parallel to the mean direc- 
tion along which the step runs (see Fig. 1) [26], while a|| is 
the atomic length unit (along (110)). We have measured 
the spatial correlation function G(y) for various angles 9 
on the Cu(5,8,90) surface at temperatures between 293 
and 328 K. For the highest temperature and larger angles 
0, the plot of G(y) showed a small influence of the time 
dependence [23, 24] of the fluctuations by the curving of 
G(y) at small y. The slope in Eq. (1) was determined for 
larger y where G(y) was a linear. The resulting values for 
the dimensionless quantity ksT/Pa^ are plotted in Fig. 
2 as gray symbols. The data points represent an aver- 
age over 100 individual steps with a total step length of 
1-2 /xm each. Included as stars are data at 9 — 0° taken 
from previous work of our group for Cu(lln) surfaces, 
n=13 and 19, for which h B T / (3{<S)a\\ = 0.019 and 0.0129, 
respectively [24, 27, 28]. 




FIG. 1: STM image of a surface region on Cu(5,8,90) where 
on the average the steps run 7.13° off the atomically dense 
(110) direction. The scan width is 24.3 nm. 



The step stiffness $ can also be determined experi- 
mentally by analyzing equilibrium shapes of islands [19]. 
Islands in equilibrium have a well-defined chemical po- 
tential which is the same all along the perimeter of an is- 
land. The chemical potential of the island is proportional 
to the product of the curvature k(9) of the island perime- 
ter and the step stiffness (3(6) (via the Gibbs-Thomson 
equation). 

For islands in equilibrium, as shown previously [19], 
the chemical potential can likewise be equated to the free 
energy per length (3o of the densely packed steps oriented 
along (110) divided by the distance of this step segment 
from the center of the island tq. Then one obtains 



P(6) = 



(3o 



t k(6) ' 



(2) 



Since the absolute number for the step free energy (3q for 
Cu(001) is known from experiment, (3q — 220 meV a^ 1 
[19, 29], Eq. (2) enables the calculation of the orientation 
dependence of fc^T '/ '(3(6)a\\ from the equilibrium island 
shape. Data were extracted from experiment by averag- 
ing over more than a thousand individual island shapes 
which were measured at temperatures between 287 and 
313 K. The temperature dependence of the equilibrium 
shapes is small in the range of 6 of interest here. Us- 
ing the mean shapes obtained for each temperature we 
have calculated the curvature as a function of the angle 
6 and furthermore, ksT/P(6)a\\ following Eq. (2). The 
temperature-averaged results are plotted in Fig. 2 as solid 
diamonds. They agree quite well with the data from the 
correlation function G(y), except for 6 = where it ap- 
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FIG. 2: Dimensionless inverse stiffness UbT j ' (3(9)a^ mea- 
sured via spatial step fluctuations (gray symbols) and via 
island equilibrium shapes (solid diamonds). Full squares rep- 
resent Eq. (6) and the dashed line is calculated from the 
temperature-dependent analytical solution of the island equi- 
librium shape in the nearest- neighbor Ising model. 



pears as if k B T / (3(0)a\\ obtained from the island equilib- 
rium shapes would be significantly higher. This is, how- 
ever, an artifact of the island shape analysis. Technically 
one cannot determine the rapidly changing curvature of 
the island equilibrium shape near the (110) direction. 

The rapid descent of k B T / (3(6)a\\ close to 6 = is a real 
effect. It can be shown that at 6 — and for moderate 
temperatures, k B T / (3(0)a\\ equals twice the concentra- 
tion of thermally excited kinks [25]: 



k B T/p{0)a\\ ^(fli/fl,,) 2 ^"/ 4 ^, 



(3) 



where e k is the kink energy. Since the unit spacing a± 
between (110) rows is the same as ay for an (001) fee 
surface, we do not distinguish between them henceforth, 
calling both a\. Eq. (3) is model independent and exact 
if z = cxp(-e k /k B T) < 1. 

In order to compare our experimental data for 
k B T/(3(6)a\ with theory, we make use of explicit ex- 
act expressions for (3(6) for the nearest-neighbor Ising 
and SOS models [16, 17], from which we obtain (3(6) as 
(3(9)+ (3" (9). While the general expression can be readily 
handled numerically, T <C e k /k B in the present problem, 
permitting an expansion in z (which is of order 6x 10~ 3 
when we specify e k below) . For the Ising model, Rottman 
and Wortis [16] found 



(3(6) ai /k B T = (e/k B T){\cos6\ + |sin6>|) 
-((I cos 6»| + |sin6>|)ln(|cos0| + |sin6>|) 



-|cos0|ln(|cos6>|) - |sin0|ln(|sin0|)) 
|cos6>| 3 + |sin6>| 3 _ 2 | 4 



cos ( 



(4) 



The energy-like first term in Eq. (4) makes no contribu- 
tion to (3(6) since it is canceled by its second derivative. 
We then obtain the following expansion for the stiffness: 



k B T/(3(6) ai = (]cos0| + |sin0|)|cos0||sin0| 
+2(|cos6>| + |sin0|) 3 (|csc0||sec6»|-l)z 2 + O(z 4 !j5) 

The first term in Eq. (5) increases gradually but mono- 
tonically with 6, while the monotonically-decreasing co- 
efficient of z 2 is sharply peaked at 6 — 0. Hence, except 
near 6 — (specifically, so long as tan# 3> z), the z 2 term 
in Eq. (5) can be neglected, leaving 



k B T 
p(6) ai 



: (tan 2 6>+| tan<9|)|cos(9| 3 = 



tan 2 6 + 1 tan 6 1 
(1 + tan 2 6)i 



(6) 



The result in Eq. (6) can be obtained from a combinatori- 
cal analysis for the number of ways of arranging the kinks 
forced by the azimuthal misorientation 6 [16, 30]. Since 
we expect k B T/(3(6)ai —^2z&s8^>0, the expansion in 
Eq. (5) evidently fails near this limit. 

The temperature independent expression for 
k B T/(3(6)ai from Eq. (6) is displayed in Fig. 2 as 
solid squares. With the help of Eq. (2) k B T / (3(6)ai can 
also be calculated for arbitrary temperatures from the 
analytical solution of the island equilibrium shape within 
the Ising model [17, 19]. The result is displayed as a 
dashed line in Fig. 2. The Ising parameter representing 
the kink energy as well as the step energy per atom was 
chosen as e fe =128 meV [24, 27]. 

While experimental data and theory agree well for 
small 6, they disagree substantially for larger 8. The 
agreement for small 6 is due to the fact that eq. (4) is 
independent of any assumption other than z <C 1. The 
reason for the deviation at larger 6 can be understood 
best by considering the special case 8 =45°. From Eq. 
(6) one obtains 



k B T/P{45°) ai = 1/V2. 



(7) 



This result can be derived stragihtforwardly by con- 
sidering the spatial correlation function G(y) for a freely 
meandering steps in the (100) direction if overhangs are 
excluded and the statistical weight factors for micro- 
scopic paths take only nearest neighbor broken bonds 
into account. As overhangs are excluded the result is 
correct to terms of order z 4 which is an excellent approx- 
imation for the temperatures considered here. The fact 
that k B T I (3(8)a\ is lower than predicted by Eq. (7) must 
therefore mean that the (100) oriented step (0=45°; 100% 
kinked) with one kink following another is energetically 
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favored over paths which deviate from this orientation. 
This is in accord with the analysis of island equilibrium 
shapes [19], which found that the ratio of line tensions 
/3(ioo) /Alio) at T = K is 1.24 rather than \p2. as pre- 
dicted in the Ising model. 

We note that diagonal, next-nearest neighbor inter- 
actions characterized by e nn cannot be the sole expla- 
nation. It is straightforward to show that the experi- 
mentally observed line tensions for the (100) and (In- 
directions cannot be fit simultaneously with the kink en- 
ergy to a model with nnn-pair-potentials. As for the 
temperature dependence, nnn-interactions were consid- 
ered, e.g., in the Akutsus' study of Si(001) [31]. To lead- 
ing order, their expansion factor cxp(— e^/fe^T) becomes 
cxp(— (efe + 2e„ n /fcsT). However, this change only af- 
fects the negligible z 2 term; there is no change to Eq. 
(6). Thus, the lower line tension of steps at larger an- 
gles 9 and, in particular, that of a (lOO)-oriented step 
(9=45°) and the lower corresponding spatial step fluc- 
tuations can therefore not be described by a nnn-pair 
potential. Rather they must be attributed to many body 
interactions. Such many body contributions might exist 
in the form of (attractive) kink- kink interactions and/or 
(negative) corner energies involved in the formation of 
kinks in a (100) -step. These interactions would keep a 
meandering step with a mean orientation along (100) 
close to the central path which is also microscopically 
a (lOO)-direction and in which kinks follow immediately 
one after another. Spatial fluctuations would therefore 
be smaller than in the Ising model. Models which might 
produce such behavior are under investigation [32]. 

We finally comment on the consequences of the strong 
angle dependence of the stiffness on the step meandering 
instability [6-13]. Under plausible assumptions, the most 
unstable, fastest growing mode A u oc (3 1 / 2 in a linear sta- 
bility analysis. Using then Eq. (3) for the calculation 
of X u not only overestimates the effective stiffness of the 
technically always slightly misoriented steps, it also in- 
troduces an unrealistic temperature dependence. A very 
recent theoretical study [33] of the meandering instabil- 
ity of Cu (001) did consider a step stiffness with angular 
dependence, but with behavior so different [34] from Eq. 
(6) to preclude useful comparison. Recognition of the re- 
markable angular dependence of the stiffness is a crucial 
ingredient for attempts to account for the meandering 
instabilities [9] and other statistical properties of vicinals 
and of single-layer islands. 
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